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Abstract 

A method based on the consistent use of the Green function formalism has been 
developed to calculate the distribution of the single-particle strength in odd nuclei 
■ with pairing. The method takes into account the quasiparticle-phonon interaction, 

ground state correlations and a "refinement" of phenomenological single-particle 
energies and pairing gap values from the quasiparticle-phonon interaction under 
consideration. The calculations for 121 Sn and 119 Sn that were performed in the 
quasiparticle£g>phonon approximation, have shown a reasonable agreement with ex- 
periment. The ground state correlations play a noticeable role and mostly improve 
the agreement with experiment or shift the results to the right direction. 
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1 Introduction 



As is well known the role of the quasiparticle-phonon interaction is essential 
for the description of excitations in odd-mass nuclei [1-3]. In magic [1] and 
semi- magic [4] nuclei it is possible to restrict ourselves to the approximation 
of the squared phonon creation amplitude g 2 in the propagators of the integral 
equations under study [4]. In other words, it is necessary to take into account 
at least the complex configurations lp®phonon or lhCgiphonon (for non-magic 
nuclei we will use the unified notation lqp®phonon). In the Green function 
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(GF) language this means that for magic nuclei it is necessary to solve the 
Dyson equation with the mass operator 



M= >°''>''°> (!) 

where the circle is the phonon creation amplitude g. 

The main difference of non-magic nuclei from magic ones is the necessity to 
take into account Cooper pairing in the nuclear ground state. This means that 
we should also consider the following simplest energy-dependent "anomalous" 
mass operators [4,5]: 



M« = >"> < n < M( 2 ) = < n < >°> ( 2 ) 

where the lines with two ingoing and outgoing arrows denote the " anomalous" 
GF and which are proportional to the gap. Pairing phonons are not 
included here because their contribution is probably small. 

It is implied in the following that the initial quantities of our problem are 
"observed" mean field described by a phenomenological potential, e.g. by the 
Woods-Saxon one, and pairing gap. The corresponding single-particle levels 
should be extracted from the observed excitation energies of non-magic nuclei 
(we will perform this procedure). The initial pairing gaps values are taken from 
experiment or solution of the BCS gap equation with the phenomenologically 
determined pp-interaction. 

So far as the nuclear pairing problem is solved within the BCS approach 
as a rule, the quasiparticle-phonon interaction is not considered explicitely 
and quantitatively on this level. If we consider explicitely the contribution of 
the graphs (2) we must also take them into account in the pairing problem 
and , therefore, to avoid double counting, change the phenomenological pp- 
interaction entering the usual BCS problem. This question will be considered 
elsewhere but here we use a simplier phenomenological procedure of "refining" 
the gap values from the quasiparticle-phonon interaction under consideration. 
It is analogous to refining the phenomenological single-particle energies in 
magic nuclei [6], [7]. 

In this work we study the role of the terms (2) and of complete taking into ac- 
count ground state correlations to describe the excitations of 119 Sn and 121 Sn. 
The important role of ground state correlations was shown long ago, e.g. for 
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the Ml excitations in 40 Ca [7] and 96 Zr [8]. For the excitations of odd non- 
magic nuclei this was considered in [9] but without any "refining" procedure 
for the case of the phenomenological mean field. The above-mentioned pro- 
cedure of "refining" the gap and single-particle energies' values will be also 
developed and realized here. 



2 Equations for one-particle Green functions in non-magic nuclei 

The general system of exact equations for the "normal" GF's G and and 
"anomalous" GF's and in a Fermi system with pairing has the form 
[10]: 

G = Go + G m - G S (1) F (2) , (3) 
F^ = G { h) ^F ( V + G^^G, 

where G and G ( h) are free GFs, i.e. the GFs of ideal gas, E, E«, ^ and 
are full irreducible self-energy parts (mass operators). Eq.( 3) should be 
supplemented by the equations for G^ and F^\ We use the symbolic form 
of equations very often. 

It is very natural to single out explicitely known components of the mass 
operators. So we write 

E(e) = E + M(e), ^ h \e) = t (h ^ + M (h \e) , (4) 
E (1) ( £ ) = E (1) + M {l \e) = A (1) + M (1) (e), 

S (2)( £ ) = £(2) + M (2) (£) = A (2) + M (2) (e) _ 

Here the first terms do not depend on the energy variable e. The quantities 
S, correspond to a mean field and E^, E^ 2 ^ describe a pairing of the 
BCS type. The quantities M, M«, M^ 2 \ (called further as M«) are 

not defined so far. We mean that they contain the quasiparicle-phonon in- 
teraction. Due to the fact that, as was mentioned before, the single-particle 
energies €\ and gaps A a are phenomenological, the quantities should give 
a contribution to E\ and A^. Therefore, to avoid the double counting of the 
quasiparticle-phonon interaction, we must "refine" the first terms of the sums 
(4) and to obtain (see below) the corresponding e x and A A from e x and A A . 
These "refined" quantities are denoted by "tilde". 

Taking Eqs. (4) into account the general system (3) can be transformed to the 
following equations (see the derivation in [4,5]): 
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G = G + GMG - F (1) M W F (2) - GM (1) F (2) - F (1) M (2) G (5) 
F (2) = p(2) + p(2) MG + G {h) M {h) F {2) - F (2) M {1) F (2) + G {h) M (2) G, 

(and the same for and The bare GF G, G w and F^, F< 2 > are the 

known GFs of Gorkov. In the A representation of single-particle wave functions 
they have the form: 



G x (e) = G x h \-e) = ^ + ^ (6) 

e — E x + i8 e + E\—i5 

F« = Ff = -^-( — i i ) 

A A 2£ A e-.E7 A + i<$ e + E x -i5 



u 2 



where u| = l-v\ = (E x +e x )/ (2E X ), E x = \J e\ + A\. The "refined" quantities 
e x and gaps A A will be determined in the next Section. 

The physical meaning of Eqs. (5) is that we have singled out explicitely the ef- 
fects of mean field and Cooper pairing in a "refined" form. In the hamiltonian 
approach, the latter corresponds to the Bogolubov's transformation but with- 
out refining from the quasiparticle-phonon interaction in the case of taking it 
into account. 

Let us define now the quantities M^. We will take them in the simplest g 2 
approximation of Eqs. (1, 2). (In [4] it was shown that the dimensionless 
quantity g 2 is a small parameter for 120 Sn, i.e. g 2 < 1 for the most collective 
low lying 2\ and 3{~ phonons.) Then the Eqs. (5) have the following graphical 
form: 



=^=^- + > n *'> %a ^- + ><? < 



.•-**•*••. ..*•**.. 

+ <^> <^ + < n <>^ 

(in our linearized version we need only two equations). This corresponds to the 
approximation of lqp + lqp®phonon. Eqs. (7) take into account the ground 
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state correlation completely because we use two terms ("forward and backward 
going graphs") both in Eq.(6) and in the formulae for all the mass operator 
Eq.(l, 2) (see [4]), and, of course, in the QRPA calculations of the phonon 
creation amplitude g. In the diagonal approximation for M xx , the secular equa- 
tion, which determines the excitation energies of an odd nucleus with pairing, 
is obtained from Eqs.(7) and has the form: 



Sa(6) = 



\ - G X M X + F [ X 1] M { X 2) FjpM{ h) + G X M { X 1] 



x 1V1 \ 

Pi 2) M x - G x h) M (2) 1 - G x h) M { x h) + Ml 
where our mass operators are given by the form 

M 1 (e) = Mi h \-e) = Y,(9nn 



=■(2) *AV 



(8) 



+ 



s,2 



M?\e) = M?>(e 



(2), 



£ _ e 2 - uj s + i'y e + E 2 + uj s ■ 
Ao , 1 1 



27 



-), 



s 2 2E 2 e - E 2 - oo s + e + E 2 + uo s - 



), 



here we have simplified the notations: index 1 = (Ai) = (ni,ji,li,mi), oj s is 
the phonon energy and g{ 2 is the phonon creation amplitude. 

The strength of the transition to the excited state A77 under consideration 
(spectroscopic factor) is given by 



[1 + q\r))(E\r) ± e Xr t ) 
Q X (E Xv ) 



where 



(9) 



(10) 



ix + M (even)x (E Xv ) A< 1,2) + M£' 2 \E Xn ) 

£\ V — — , L±\r] — 



1 + Qx v 



'1 + Qxv) 



(11) 



M( dd)\(E Xv ) 



E 



At; 



9 A (e) = (e - e x - M x (e))(e + e x + M [ x h \e)) - (A A + M x ' (e)) 



:i2) 
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M( even ) and Mf^m are the even and odd terms of the mass operator M [11]. 
The denominator in Eq. (9) is the energy- derivative. Eqs.(8-12) determine 
the characteristics of odd nuclei with pairing in our approximation. 



3 The refinement of phenomenological single-particle energies and 
gaps 



The initial general system (3) can also be transformed to another form (see 
the derivation in [4]). Let us introduce the GF 



G = G + G (t + M)G = G + GMG (13) 

where G determines "refined" quasiparticle energies e\ (we mean Landau's 
quasiparticles here). Then the system (3) can be written as one equation for 
G: 



G — G — GTS^G^TP^G. (14) 

In Eq.(14) we use the approximation which is diagonal in the single-particle 
index A. Let us represent the mass operator M as a sum of its odd M( ocM ) 
and even M( even ) terms [11] and determine the energies E\ of the observ- 
able quasiparticle levels as dominant solutions of Eq.(14). Then we obtain the 
general formulae which connect the observable {^a, Aa} and refined {e\, Aa} 
quantities: 



even)X 

£X = 1+^) ' (15) 

. .(1,2) _ A^+Mj^jE,) 



Aa = A 



i + qx(E x ) 



where E x = yJel + A\, q x = -M ( oc id)\(E\) / E\. Thus, if the phenomenological 
quantities E\ and Aa are known we can find the bare ones, which enter Eqs.( 5), 
( 6), ( 7), from the solution of the nonlinear relations (15). 



4 Calculations of single-particle strength in ng Sn and 121 Sn. 



At first we developed and realized the procedure to extract the phenomeno- 
logical E\ and A a from the observed excited quasiparticle levels E\. We used 
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Table 1 



Refined neutrons single-particle energies i\ and gaps A a for 120 Sn. 



A 


£\, 
Mev 




A A 
MeV 


Aa, 
21 phon. 


7A % 


21 phon. 


3 phon. 


21 phon. 


3 phon. 


2d5/2 


-2.67 


-3.69 


-3.41 


1.35 


-0.23 


117 


85 


lg7/2 


-1.36 


-3.27 


-2.54 


1.58 


0.76 


52 


33 


2d3/2 


-0.09 


-0.17 


-0.24 


1.36 


1.09 


20 


18 


3sl/2 


0.13 


0.37 


0.20 


1.27 


0.89 


26 


16 


lhll/2 


0.42 


0.86 


0.91 


1.62 


1.41 


13 


9 



the formula E\ = A + E^ x , where A is the odd-even difference and E* x are 
the excitation energies of n9 Sn and 121 Sn, and determined the quantities e\ 
and Aa from an iteration procedure. This procedure was organized in such a 
way that the E\ values could be reproduced by the formula E\ = \Je\ + A 2 . 
The E\ and A A obtained are given in Table 1 (for detailes see [4]). 

Further, it is necessary to perform the "refining" procedure, i.e. to find E\ and 
Aa from the known E\ and Aa. To do it the non-linear Eqs.(15) with our g 2 
choice of have been solved. The results are given in Table 1 where 7a and 
7 quantities were defined as follows 



7A A a ' 7 Ea(2ja + 1) ' 1 J 

These values give the contribution of the quasiparticle-phonon pairing mech- 
anism caused only by the retarded pp-interaction, i.e. by the con- 
tribution in Eq.(15). In Table 1 the calculations using 3 low-lying collective 
2+, 3i , Ai phonons and 21 phonons are given for 5 levels near to the Fermi 
surface. The phonons were calculated within the standard theory of finite 
Fermi systems with phenomenological E\ and Aa which is permissible in our 
g 2 approximation because we omitted the g A terms. We obtained that the 7a 
values depend rather strongly on A and the mean value 7 = 32%, which was 
calculated for all the 8 levels under consideration. 

At last, using the E\ and Aa values as the initial data we solved Eq.(8) and 
calculated the spectroscopic factors, Eqs. ( 9, 12). The results and their com- 
parison with the experiment [12], [13] are given in Table 2 for all the 5 levels 
below and above the Fermi surface. We obtained a reasonable agreement with 
experiment except for the 2d5/2 and 1^7/2 levels in 121 Sn where strength is 
very fragmented and small, i.e. it was not observed very reliably. 

In Figs. 1,2 and Table 2 the results "GSC-" obtained without taking into 
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Table 2 

Energies and spectroscopic factors of neutron states in 119 Sn (upper line) and 121 Sn 
(low line)*. 



A 




Sx 


vl 






V.. 1 kJ v.. 


OCS— 

V.I V. l 


Exp 






u 2 , 

A 


2o & 


7 1 n 


7 


6 Q1 

U. i/ -L 


9 49 
(4.3-10.2) 


9 1 S 


S ^6 


1 



















lo9/2 


5.41 


5.90 


5.40 


0.53 
(3 9-6 51 


0.43 


0.43 


S3 1 
















S3 


2d5/2 


i i n 






n 4^ 


n 4fi 


n 3 


0.95 




1 14 


1 41 


1 66 


n 1 1 

u. ± ± 

(0 0-2 01 


n D4 


n D9 


0.05 


lo7/2 


n 7Q 


1 81 


95 


75-0 6 

\J • i <J j u. u 


66 

U. \J\J 

(0 0-1 0) 

I U.U O.U 1 


77 


0.83 




2.85 


1.23 


0.94 


0.15 
CO 0-3 51 


0.04 


0.16 


0.17 


2ci3/2 


0.02 


0.01 


0.02 


0.4;0.45 


0.4 


0.48 


0.52 







w 


w 


0.44;0.65 


0.35 


0.43 


0.48 


3al/2 


n 


~ 


~ 


96-0 39 


36 


0.42 


0.45 




0.06 


~ o 


~ o 


0.3 


0.43 


0.51 


0.55 


l/ill/2 


0.09 


0.18 


0.15 


0.29 


0.26 


0.34 


0.37 




u.uu 


n 1 7 

U. J. I 


~ o. 


n 4Q 


n 4s 

U.t:0 


n ^6 


0.63 


2 f7/2 














~ 




2.83 


3.11 


2.7 


0.35 
(0.0-4.0) 


0.39 


0.59 


S3 1 


3p3/2 














S3 




3.73 


4.97 


4.9 


0.54 
(0.0-5.0) 


0.42 


0.4 


1 


l/i9/2 

















+H13/2( & ) 


7.5. 


7.47 


6.61 


w 24 


19.92 


22.23 


1 



*For the states without strong dominancy the energy intervals of summation are 
given. The energies E x correspond to the dominant peak position or to the mean 
energy. 

& The quantities S\ * (2j + 1) are given, g 



SiE^MeV- 1 




EfP = hAlMeV 



UE x (MeV) 



Fig. 1. Strength distributions for the 1^9/2 neutron state in 119 Sn. The solid line 
gives the results with taking into account all the ground state correlations, the 
dashed line gives the results without them. The averaging parameter is l.OMeV. 
The experimental data (histogram) are taken from [12]. 
S(E x )MeV~ 1 

0.7' 



0.6 
0.5 
0.4 
0.3 _ 
0.2 _ 
0.1 _ 




Ef p = 2.83MeV 




n 6 » 10 12 \±E x (MeV ) 

Fig. 2. The same as in Fig.l but for the 2/7/2 neutron state in 121 Sn The experi- 
mental data (histogram) taken from [13]. 

account all the ground state correlations (except for the QRPA ones in our 
phonons) are also given. We see that the difference is rather noticeable. The 
inclusion of GSC mostly gives a correct trend, both for the levels near to the 
Fermi surface and for those beyond this surface, and improves the agreement 
with experiment. One can see from Fig. 1,2 that it is more manifested for the 
"differential" results than for the integral ones presented in Table 2. 
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We also investigated the role of the terms M ^ , Eq. ( 2) of the anomalous 
mass operators both for the ground state ("refining" the gap) and for the 
description of excited states. It turned out that they give small contributions 
only for states corresponding to the single-particle levels which are very far 
from the Fermi surface. For the rest of the states considered, they have a 
contribution and improve the agreement with the experimental data rather 
often. 

In conclusion, in order to calculate the distribution of single-particle strength 
in odd nuclei with pairing we formulated the method which takes consistently 
into account all ground state correlations (within the approximation used) and 
new terms, i.e. Eq.(2), of the "anomalous" mass operators which are specific 
for non-magic nuclei. For the first time the method includes also the refinement 
of the phenomenological pairing gap values from the quasiparticle-phonon in- 
teraction under study. A more exact than usual receipt of determing phe- 
nomenological single- particle energies for non- magic nuclei was realised here. 
The first calculation for 119 Sn and 121 Sn performed in the lqp + lqp®phonon 
approximation showed a reasonable enough agreement with the available ex- 
perimental data. 

The mean value of the contribution of the quasiparticle-phonon pairing mech- 
anism caused only by the retarded pp-interaction was obtained for the first 
time, it is about 32% for 120 Sn. A noticeable numerical role of ground state 
correlations, although not so dramatic as in [9], was obtained. 

S.K. thanks Profs. P. von Brentano, U.Kneissl and Dr. N.Pietralla for very 
useful discussions. 
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